VOLUME SEPARATE No. 223 


PROCEEDINGS 


SOCIETY 


ENGINEERS 
JULY, 1953 


APPLICATIONS THE RELAXATION 
TECHNIQUE FLUID MECHANICS 


John McNown, En-Yun Hsu, 
and Chia-Shun Yih 


Presented 
Miami Beach Convention 
June 16-19, 1953 


ENGINEERING MECHANICS DIVISION 
{Discussion open until November 1953} 


Printed the United States America 


Headquarters the Society 
39th 
New York 18, 


PRICE $0.50 PER COPY 


AMERICAN 
; SOCIETY OF 
| 
ENGINEERS 
POUNDED 
1062 


THIS PAPER 


--represents effort the Society deliver 
technical data direct from the author the 
reader with the greatest possible speed. 
this end, has had none the usual editing 
required more formal publication procedures. 


Readers are invited submit discussion apply- 
ing current papers. For this paper the final 
closing dead line appears the front cover. 


Those who are planning papers discussions 
for “Proceedings” will expedite Division and 
Committee action measurably first studying 
the printed directions for the preparation 
ASCE technical papers. For free copies these 
directions—describing style, content, and for- 
mat—address the Manager, Technical Publica- 
tions, ASCE. 


Reprints from this publication may made 
condition that the full title paper, name 
author, page reference, and date publication 
the Society are given. 


The Society not responsible for any statement 
made opinion expressed its publications. 


This paper was published 1745 State Street, 
Ann Arbor, Mich., the Society 
Civil Engineers. Editorial and General Offices 
are West Thirty-ninthStreet, New York 18, 
N.Y. 


APPLICATIONS THE RELAXATION TECHNIQUE 
FLUID MECHANICS 


Relaxation, numerical method integration, can 
used obtain solutions for wide variety problems which 
are not solvable standard employing this meth- 
od, network values assumed for the function sought, 
and then systematically corrected the errors are "relax- 
ed", Based the calculus finite differences, this method 
was brought fruition Ingland Southwell and his assoc- 
lems elasticity, heat transfer, and fluid mechanics during 
recent years. Although only particular solutions can 
tained such methods, and although the computations are time 
consuming, relaxation methods are effective cases for 
general solution cannot obtained direct methods. 


Partial differential equations such the Laplacian, 
the Poisson, and the biharmonic can integrated using the 
relaxation procedure. These equations are often encountered 
engineering, and general solutions have been 
tained for simple boundary conditions. However, general solu- 
tions are usually unobtainable the boundary conditions are 
complex. such cases that numerical solutions are 
particularly valuable because their application complex 
boundary condition can fulfilled almost easily 
ple one. this discussion presented only the method 
solving the Laplace equation, the most frequently encountered 
the analysis fluid flow. Although this the simplest 
the several equations, the manner solving nonethe- 
less indicative the general method, Applications this 
differential equation are potential flows such 
seepage, flow through boundary contractions, 
slots orifices, certain types wave motion, and flow over 
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weirs. for and axisymmetric flow 
terms feasible obtain solutions. 


Resolution any significant problem means the 
relaxation process will require least many hours and per- 
haps long several months, yet the individual steps 
the process are few and easily These steps are pre- 
sented some detail, herein, from the initial assumed 
values, through the elimination residual and the 
advances successively finer nets, the eventual satisfac- 
tion both the differential equation and the assigned bound- 
ary conditions. number applications are discussed the 
last section, the analyses having been conducted the Iowa 
Institute Hydraulic Research graduate students and mem- 
bers the staff. 


Much the described herein was supported the 
Office Naval Research under contract Among 
those who assisted the tine-consuming computations are A.H. 
Abul-Fetouh, A. LeClerc, Ze C. Soong, WwW. Hsieh, s. Ince, Ss. C. 
Liu, and Sarpkaya, either now formerly graduate students 
and research assistants the State University Iowa. Dr. 
Hunter Rouse, Director the Institute, suggested directed 
several the studies. 


Essential Elements Finite Difference Theory 
The goal applying the relaxation method the 


tainment network values function which satisfies 
given partial differential equation each point dis- 


Fig. Pattern function values; 
(a) regular, and (b) irregular 


crete network covering the region, and which also satisfies 
given set boundary conditions. number points must 
large enough (or the mesh fine enough) that the assumption 
linear variation the function between adjacent points 
justifiable. For such network, the theory finite dif- 
ferences can successfully used. 
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algebraic expression for the Laplacian partic- 
function, say can obtained quite simply. With 
Fig. the approximate value 


sought for the central point terms the several values 
the immediate vicinity, the mesh being square and the 
lines spaced distance apart. the coordinates the 
central point are designated and the points 


respectively, The second derivative can now 
obtained similar manner: 


the same way, one obtains expression for the second par 
tial with respect 

Finally, the equivalent form for Eq. (1) expressed finite- 
difference form is: 


Although Eq. (2) approximation, equality ap- 
proached the spacing For given set 
values, either assumed, can evaluated 
means Eq. (2), and its difference zero (or other 
signed value for the Poisson equation) noted. The successive 
elimination these undesired residues discussed the 
following section. 


Corresponding relationships can also derived for 
flow, but additional complication 
The appropriate form for the differential equation the 
locity potential follows: 


(3) 


which measured parallel the axis and 
right angles it. Although both the velocity potential 
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and the stream function for two-dimensional flow satisfy 
the sane equation, the Stokes stream function, for axisymmetric 
flow, satisfies equation similar Eq. (3), but differing 
the sign the first-order term: 

(4) 


case (Fig. la), 


this value substituted into Eq. (3), together with the 
finite-difference the second partial deriva- 
tives, 


the corresponding form for 


the axis symmetry and the point 


the four neighboring points are not equidistant 
from given point, for point near irregular boundary, 
the forms Eqs. (2), (5), and (6) must 
shown Fig. lb, the lengths two the legs the star 
may less than these lengths are designated 

and the appropriate form for the La- 
piacian operator for two-dimensional flow 


and for axisymmetric flow, 


262 gq, | g, 
For some problems preferable use reverse 
method which the dependent and independent variables are 
interchanged. two-dimensional potential flow, and 
are considered the independent variables, and 
satisfy the Riemann-Cauchy equations 
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from which follow the inverse Laplace equations: 


2 2? 


These equations can, course, expressed terms fi- 
nite differences the aforanentioned manner, 


Because the region interest fluid flow can al- 
ways bounded pair potential lines (on which 
constant) and pair stream lines (on which constant), 
irregular stars will encountered the plane 
chosen for computation. This advantage sometimes compensates 
for the drawback that the values and along sme 
are unknown (though related) functions 
illustrated later, the use and independent 
variables particularly advantageous if, prob- 
the boundary lines are not fixed priori, but 
determined trial from double boundary condition. 


Technique 


The relaxation process can described method 
addition the special techniques for the satisfying 
various boundary conditions, two principal procedures are 
(a) the reduction the residues, and (b) the sub- 
iivision the network wherever the linear 
iation found inaccurate. 


The region flow first drawn large scale and 
then subdivided into and approximate values the 
function selected each intersection. For example, the 
stream function can from sketched stream lines 
once arbitrary values this instance) have 
been assigned for the center line and the boundary. The nu- 
merical values can selected will, because the results 
are presented dimensionlessly terms the char- 
acteristics either upstream downstream from the transition, 


Using Eqs. (2) and (7), one can compute from the as- 
sumed values the function the value the residue for each 
the intersections. For example, near the point marked 
Fig. the values could shown Fig. and the 
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idue recordea the right the intersection the lines. 

other values are computed. order reduce 


dimensional transition 


the necessary adjust repeatedly the many 
values 


For regular star, one for which all four legs are 


the corresponding residue 
apparent fran (2), 
(6). Thus for Fig. 
residue reduced from 
ring (2). the same 
adjustment, the residues for 
each the neighboring points 
are also decreased 
cause becomes the 
computation the residue 
the point The process 
cordance with the simple star 
Calculation patterns shown the left 
initial residues Fig. for the several possi- 
ble cases. Because each 
change affects four other directly and still 
others indirectly, the various values will need correct- 
many Large residues are first reduced without at- 
tempting make them extremely small, refinement being sought 
only after the large discrepancies have been eliminated. 
less than hour computation one can obtain familiarity 
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Two-dimensional 


Fig. Patterns for relaxation 


with the simple process and feeling for various short cuts 
such and regional block relaxation for 


process (at points 


which simple formulas can derived. Fig. shown 
stage calculation following single operations each 
points (in alphabetical order). The original 
were obtained from Fig. the values having been multiplied 


223-7 


= 
(b) Axisymmetric pottern 
-18 
38 


avoid the use decinals. 


Reduction the residues negligible amount 
coarse net such shown Fig. usually not sufficient 
because errors are likely made the many numerical 
processes and because the net too coarse for the assumption 
linearity hold, The errors can found recamputing 
all the residues using the corrected values the function, 
difference between the recomputed value and that obtained from 
the successive alteration indicating error calculation 
bookkeeping. Thus for point Fig. the calcula- 
tion 


based the revised values constitutes check several 
the preceding calculaticns. 


Such errors once found can usually eliminated rath- 
rapidly and hence give particular trouble unless resi- 
dues are checked too they occur dis- 
tance from the boundary flow pattern, 
corrections can made directly accordance with 
termined pattern (Fig. 6). the preparation Fig. 
was assumed that residue (or error) 100 was found for the 


Fig. Reduction single residue 


point and that those for all other points the vicinity 
were negligibly small. reducing zero, 128 
the neighboring values the function were found 
other than 100, proportional values can used, course. 
The correction any point determined directly terms 
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the the points; point three squares 
the left and two squares below the point which the error 
was made, the correction indicated Fig. would 
the anomalous 


Because the initial grid deliberately made coarse 
for simplicity, the variation some regions 
ably not linear. The existence this undesirable condition 
can easily ascertained once the initizl has been 
satisfactorily completed. given region value can 
obtained for the center square averaging the four 
this done for each four squares with vertex, 
the four new values can then averaged second val- 
for the common this say differs 
significantly smaller spacing should used. 
The necessary step, known advance firer net, need 
taken cnly the region for which evident 
from such computetions. successively finer nets are used 
only required, the amount work kept minimun, 
Fortunately, the check for and the first step 
advance are the same. the work progresses, 
able add one more zeros each the numbers, thereby 
further reducticn the residues without resorting 
amount experience required insight into the 
proper for balancing the merits additional 
one. Naturally, the relaxation technique that has been de- 
scribed applies alsc and are used independent var- 


The accurecy given result can only judged 
computing velocity pressure distributions along boundary 
cther arbitrarily selected line and determining whether 
they are both systematic and comparatively insensitive 
ther For this reason desirable make such 
computations various stages. may known that small 
region non-linearity remains, but further advance does 
not aiter appreciably the calculated distribution nearby, there 


Conditicns 


Little mention has been made the method 
fying given set boundary conditions because the procedure 
actually separate part the problem, and because 
iety techniques required. the transition Fig. 
difficulty exists. The values along the solid bouni- 
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ary and center line are fixed, and the net simply, labo- 
riously, extended the right and left the condition 
constant velocity ky) satisfactorily approached. 

For other problems, however, special procedures are required. 
Among such conditions are those specifying (a) symmetry the 
boundary, (b) free streamlines with without gravitational 
effects, (c) the free surface seepage problem, and (d) 
the special condition for the inverse method. some cases, 


Fig. Computation Fig. Indication 


the velocity the boundary must evaluated before the 
boundary condition can 


Along solid boundary, the function variable 
and aconstant. used, the velocity can com- 
puted directly for the midpoint each interval along the 
boundary from the simple formula 


= 


which represents distance along streamline and 
and are values along the boundary shown 
Fig. lb. used the normal gradient 
mated from two values, one the boundary and the other near 
the boundary. such case (see Fig. 7), the calculation 


oy 


yields approximation the the velocity 
point which distance A,a/2 away from the boundary. 
Even though the assumption local linearity approaches the 
actual conditions, better results are obtained using 
finite-difference formula for the derivative the boundary 
based three W-values illustrated Fig. 
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Inclusion more than three terms was found unnecessary, 
the increased complexity the computation being 
panied significant increase accuracy. Curvilinear 
trapolation this kind might sean unwarranted method 
which based the assumption locally linear variation, 
but the resulting increased accuracy can obviate advance 
finer net. 


the event that flow pattern symmetrical with 
respect center only one the two halves need 
determined, the line symmetry being Fig. 
automatically fixed for solid and resi- 
dues which would transferred points such boundary 
are ignored. If, however, the velocity 
used, the boundary values are not known but are subject the 
that lines constant must intersect the line 
symmetry right angles satisfy this re- 
only necessary note that any residue coming 
the center line from point the neighboring 
line (designated Fig. would also come from the 
point svmmetrically placed the other side the 
center, and therefore All other steps 
the process remain unchanged, points line and the 
cer.ter line teing relaxed the usuel 


the analysis the patterns flow occurring for 
verious types efflux, the free streamlines which are 
countered necessitate special measures. Along free 
line the pressure constant; thus, for irrotational flow, 
the velocity either constant, gravitational effects are 
ignored, proportional the square root the 
distance the line tctal head, gravitational effects 
are included. Initially the location the free streamline 
not known, ana must assumed. The 
ponding flow pattern next through use the re- 
laxation process, and the velocities the free stream- 
lines are the obtained will not, 
general, satisfy the specified boundary condition for the 
velocity, the assumed location the free streanline must 
changed and the process Because used for the 
direct this type problem, the aforementioned 
extremely 


the analysis unconfined seepage porous 
medium, somewhat different type free 
countered. the piezometric -head (velocity 
generally used, two conditions are required: (a) every 
point along the free streamline the pressure must zero, and 
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(b) the normal canponent velocity must zero. 
ly, after trial curve for the streamline bounding the flow 
hes been drawn, values can assigned for which automatic- 
ally satisfy the condition that that is, that 
equal the After the corresponding complete 
work has been obtained, the second boundary condition can 
checked, and revised bounding line drawn required. Once 


Fig. Schematic diagram for reverse method 


again, this process can ropeated, entirely 
until acceptably correct result obtained. 


check the imposed boundary condition for 
ating the velocity along curved irregular boundary, 
Fige determined from (11) for the point 
the bounding streamline. The value can then 
determined from the expression 


which the angle between the bounding streamline and 
the x-axis. For axisymmetric flow, the corresponding expres- 
sion 


(No sign convention need defined because the direction 
the velocity never question.) 


The treatment the boundary condition and 
are used independent can best explained 
means example. the surface jet issuing from 
two-dimensional nozzle, the pressure and magnitude the 
locity are constants. This condition will determine the shape 
the jet. With reference Fig. sufficient 
consider the region ABCDE, where and are comparatively 
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far away from the contraction. Along CD, can assumed 
constant, say zero. The corresponding value for can 
along ED, one has along ABC, being the discharge 
per unit width the The value along can 
taken sufficiently large constant that the flow 
downstream from the contraction essentially uniform, but 
should remembered that far the corresponding value 
(denoted has not been Aside from the 
somewhat complicated boundary conditions for and 
ABC, the conditions the boundaries are: 


Along would possible translate the 


being measured direction normal AB) into one 
volving the unknowns, but the result extreme- 
unwieldy. Instead, one demands the satisfaction one 
the equations (Eq. (9)): 


and requires that 
and 


equivalent 
f(x) 
which describes the boundary AB. 
The conditions along can similarly formulated, 


Qne again demands the satisfaction one the relationships 
from Eq. (9) along But instead Eqs. (12), one 


requires 


which the ratio the ultimate width the jet 
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With reference Fig. 10, the procedure then 
follows. Assuming along AC, one obtains 
obtained integrating (14) numerically, being 
sumed such that the calculated will equal 


x-plane y-plane 


Fig. 10. Work sheets for the ‘evaluation 
and 


With the value given, the boundary con- 
dition for determined from the linear rela- 


tion between and along CD. 


After the relaxation procedure for campleted, 
the result checked against one Eqs. (9). this equa- 
tion not satisfied, corrected accordingly. The 
computed either one relaxation process, directly 
using Eqs. (9). 


approximately -ad/q and approximately 
approximately equal The final result computa- 
tion for the two-dimensional nozzle with the boundary geometry 
given Fig. presented the following section. 


Once the various steps the relaxation process, 
outlined above, have been completed, any required results can 
obtained: the shape free streamline established 
the process, the pressure and velocity distributions can 
determined throughout space with reference the pressure and 
velocity region uniform flow, comparative quantities 
seepage can evaluated, and coefficients discharge can 
computed. The one remaining question the degree ac- 
curacy obtained difficult one answer. 
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There direct way assess the accuracy other 
than make determinations of, for example, velocity 
from each the various nets prepared. the change 
this result following particular refinement 
tively small one can conclude that the result satisfactory. 


Fig. 11. curved inlet (S. Ince) 


Various ways resolving partially this troublesome matter 
will become evident each significant fact 
that two regions are entirely equivalent. one region 
several subdivisions may required, another none. 
thermore, even sizeable errors the pattern one point may 
have very small effect another. the entire process 
time consuming, special care should taken assess each 
stage the calculation the reliability the results. 


Applications 


The importance the relaxation method best indi- 
cated its diversity application which evident from 
the following series examples for flow through boundary 
transitions, flow with free surface, and seepage. The 
terms are axisymmetric. the details 
the procedure for the interior are similar all 
problems, only the satisfying the boundary conditions 
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the problems differ significantly. 


The results calculation the flow through 
two-dimensional inlet transition between reservoir and 
conduit are shown Fig. 11; this instance the radius 
rounding equal half the conduit width. The distribution 
pressure along the boundary indicated, representa- 
tive streamlines are shown together with the radial lines 
which they approach asymptotically. The pattern the flow 
some distance upstream the inlet identical that for 
flow toward sink the intersection the projection 
the vertical boundary and the axis symmetry. The result 
has been found conform closely with that obtained from 
different type analysis and with that fram experiment 


second example internal boundary transition 
axisymmetric flow. study series head forms [3], 


1.0 


-0.6 


Fig. 12. Pressure distribution for axisymmetric flow 
past modified ellipsoidal head form (W. Hsieh) 


the pressure distribution was determined for modified el- 
lipsoidal form indicated Fig. 12. comparison made 
the results obtained from the Institute water tunnel with 
those calculated numerically. The correspondence good from 
the point stagnation point somewhat beyond the region 
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0.4 


minimum pressure. The dimensionless value for the minimum 
pressure -0.57, which can compared with value 
for hemispherical nose form. The significance this 
parison that body with the modified nose form could move 
higher without producing cavitation than could 
one with the hemispherical form. The effects the 
layer growth this regicn are evidently very small. 


Applications the relaxation technique the 
mination flow patterns involving free surface include 
problems efflux, overflow, jet impingement, and steady 
ity flow. these characterized the existence 
boundary along which the pressure constant. These have 
been the subjects several calculations which serve 
ther illustration the 


From extensive calculations, Rouse and Abul-Fetouh [4] 
have presented detailed information for efflux through ori- 
fice symmetrically placed the end circular conduit. 
The shape the free jet, the pressure distribution along the 
solid boundary, and the ccefficient contraction were pre- 
sented for various ratios the diameter the orifice 
the diameter the conduit. comparison the results with 
those for comparable flows revealed remarkably 
close correspondence for the values the contraction coeffi- 
cient. Furthermore, the coefficients available tron various 
experiments were found follow essentially the same trend. 


Overflow problens are inaccessible ordinary 
ods hydrodynamics because the paramount importance 
gravity. That is, surfaces constant pressure are not, for 
these flows, also surfaces constant velocity. even 
the classical method the theory not 
plicable. Southwell has presented calculation for the 
free Also, Citrini recently published the 
sults calculation for circular weir which are appli- 
cable the design spillway; obtained 
the pattern for the case which (horizontal) radius 
curvature the weir ten times the total head the weir. 


The results refined calculation for 
sional flow over very high weir, conducted the Iowa 
stitute Hydraulic Research, are shown Fig. 13. 
again, the flow great distance the left could taken 
that resulting from sink the intersection the plane 
the weir with the surface total head. Far the right the 
streamlines approach trajectory free fall, and the inter- 
nal pressure approaches zero. The significant region near 
the crest the weir shown three different scales 
the various parts the figure. The coefficient discharge 
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the conventional equation 


3/2 
(15) 


was found 0,408. This value somewhat more, 
should be, than the value 0.388 found Citrini for the 
curved weir. Both free surfaces and representative internal 
streamlines are shown. the shape the former 


= q 
v¥2qg 
4 


pressure-heo 
4/4 
Fig. Pattern flow over 


very high weir Liu) 


cant the design spillway, coordinates the two bound- 
ing curves are given the table their ratio 
the 


Perhaps the longest all the calculations described 
herein, performed the second author this paper part 
his doctoral dissertation, was that for axisymnetric flow 
internal boundary which was partly fixed and partly 
steady-state cavity. Once again, this calculation duplicated 
part the experiments described reference being for 
flow around hemispherical head form mounted shaft 
equal The cavitation index 
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Po- 


for this case was 0.2, sufficiently iow that large 
vapor cavity formed. The pressure distribution observed 
the experiments was imposed internal boundary condition, 


Coordinates Free Surfaces for 


Flow over High weir 


Upper 


Nappe 8x/H 


0.991 
0.989 1/32 
0.988 1/16 
0.985 
0.983 
0.973 
0.969 
0.965 
0.959 
0.951 
0.943 
0.932 
0.921 
0.903 
0.880 


and the corresponding shape the vapor cavity was determined 
successive 


Fig. are shown (a) the observed and 
cavity profiles, (b) the pressure distribution, and (c) the 
pattern used for the calculation. details the cal- 
culation are available thesis form [6], including the 
values each the intersections Fig. The condi- 
tions far upstream far downstream were, course, those for 
parallel means artifice, the radial extent 
the net was restricted ten times the radius the body. 
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Upper Lower 
0.855 
0.920 
= 0.031 
0,947 
0.068 
0.731 0.113 
0.735 0,101 
0.965 
0,609 0.019 
0.531 =-0.958 
0.422 -0.135 
0.340 
0.230 
0.114 -0.432 
-0.910 


axis. This was the topic investigation the lowa In- 
stitute LeClerc [7]. this study, the form the free 
surface was determined electrical analogy, only the internal 
flow pattern being evaluated relaxation. The velocity was 
assumed sufficiently high that the effect gravity 
well those viscosity could neglected. The shape 
the jet, representative streamlines, the variation the pres- 
sure along the axis and along the plate, and the liniting curve 
for the free surface are shown Fig. 15. 


characteristic application the inverse method 
shown Fig. 16. The pattern efflux from curved two- 


Fig. 15. Impingement jet flat 
plate (A. LeClerc) 


dimensional nozzle that which was also used explaining 
the techniques. arbitrary parabolic curve the end 
parallel approach reduces the normal section the region 
efflux one third that the approach. contraction 
coefficient was found 0.75. point which worth 
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14. Flow with cavitation for hemispherical 
head fom (a) cavitation pocket, (b) pressure 
distribution, and (c) relaxation pattern 
(E. Hsu) 


The values some distance (radially) away fran the body 
were computed for axisymmetric half body equivalent size 
and location. was then necessary extend the net only 
until its values coincided with those calculated for the half 
body. The resultant saving time was considerable one. 


third type free-surface flow typified the 
impingement circular jet plate normal the jet 
223-21 


RELAXATION 
RELAXATION 
EXPERIMENTS 
— = — — 


peating that regular stars are used throughout this calcu- 
lation even though both the solid and the free boundary are 
curved, 


Problems seepage third group which 
the methods relaxation are applicable. well known, 


Fig. 16. from curved 
nozzle (T. Sarpkaya) 


consequence the Darcy law the applicability potential 
theory problems laminar seepage. The pressure the 
piezometric head acquires the role the velocity potential, 
and the Laplace hence applicable. examples 
the application the relaxation procedure various types 
seepage flows have already been published. Results 
ply bounded structures have been presented one the 
authors [8] for two-dimensional horizontal seepage through 
monolith, and McNamee for various patterns two- 
dimensional seepage into excavation which partially pro- 
tected impervious sheet piling. 


references [1] and [8] are given applications 
the more complex case seepage with free 
version the latter being presented Fig. 17. These 
are cases for which, once again, the streamline 
initially unknown, and can treated only trial and 
The two boundary conditions for the free surface are: (a) the 
direction the resultant velocity must coincide with that 
the free streamline, and (b) the free surface must one 
constant pressure (i.e., the piezometric head, 
equal the The second these was assumed 
for trial curve and relaxation the interior was 
Adjustment the free surface was then made the basis 
the shortcomings with regard the first condition. the 

was used throughout, the boundary condition for 
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the horizontal bottom line was that orthogonality (or 
and was fulfilled the application double res- 
idues already described. The value constant 
the upstream face, but equal the distance below the 
section the free surface and the downstream face the 
latter. 


Although many problems seepage two dimensions 
can obtained more satisfactorily means classical 
thods, only such method relaxation can used for 
plex curved boundaries. addition, axisymmetric 


Fig. 17. Seepage through vertical 
wall (C. Yih) 


such flow ground water well the axisymmetric 
counterpart problem seepage into excavation, 
are beyond the scope ordinary 


Resume 


the presentation both the relaxation method and 
examples which utilized, the utility and flexibili- 
the process have been amply demonstrated. 
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0.43h 
yh 


the technique can quickly mastered, containing does 
only two basic parts. The mechanism the internal reduction 
residues simple and extremely repetitive. The various 
boundary conditions are either assigned initially met in- 
directly part the process. advantage the method 
that absolute accuracy not requirement; occasional 
error once made will caught and need not introduce ex- 
cessive loss time. Although the method simple, 
not completely mechanical, that the exercise judgment 
and the introduction suitable modifications given 
tine can greatly reduce the time required. 


The method restricted the resolution problems 
essentially irrotational flows for flows for which the 
Reynolds number quite small. will not and should not 
considered replacement for the more general analytical 
solutions characterizing classical Although 
still restricted flows which are essentially two 
axially symmetric, its utility lies the fact that 
within these limitations the relaxation process can used 
regardless the complexity the boundary condition. 
small change boundary form that might preclude the solving 
otherwise classical problem standard methods presents 
difficulty its resolution means numerical methods. 
the engineering profession the process 
ing many new applications for analyses potential flows, the 
relaxation technique recommended powerful adjunct 
standard methods. The flow patterns, presented herein for the 
first time with only two exceptions, are both useful 
selves and representative the variety problems which can 
solved relaxation. 
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